cc// ' 


NASA Contractor Report leeooi 


NASA-CR-166001 

19820023830 


THE CRACK PROBLEM FOR A NONHOMOGENEOUS PLANE 


F. Delale and F. Erdogan 


LEHIGH UNIVERSITY 
Bethlehem, Pennsylvania 18015 


Grant NGR 39-007-011 
September 1982 


LIBRARY COPY 

OCTO 1982 


LftNGLEY RESEARCH CENTER 
LIBRARY, NASA 
H.AWPTCN, VIRGINIA 


IW\SA 

National Aeronautics and 
Space Administration 

Langley Research Center 

Hampton, Virginia 23665 



NF01898 



THE CRACK PROBLEM FOR A NONHOMOGENEOUS PLANE* 

by 

F. Del ale** and F. Erdogan 
Lehigh University, Bethlehem, PA. 18015 


Abstract 

In this paper the plane elasticity problem for a nonhomogeneous medium 
containing a crack is considered. It is assumed that the Poisson's ratio of 
the medium is constant and the Young's modulus E varies exponentially with 
the coordinate parallel to the crack. First the half plane problem is for- 
mulated and the solution is given for arbitrary tractions along the boundary. 
Then the integral equation for the crack problem is derived. It is shown 
that the integral equation having the derivative of the crack surface dis- 
placement as the density function has a simple Cauchy type kernel. Hence, 
its solution and the stresses around the crack tips have the conventional 
square-root singularity. The solution is given for various loading condi- 
tions. The results show that the effect of the Poisson's ratio and conse- 
quently that of the thickness constraint on the stress intensity factors are 
rather negligible. On the other hand, the results are highly affected by 
the parameter 6 describing the material nonhomogeneity in E(x) = Eoexp(3x). 


1. Introduction 

In practical applications the material nonhomogeneity becomes an impor- 
tant factor to be considered particularly in two classes of problems. The 
first is a group of problems in geophysics in which, because of the size 
of the medium, the spatial variation of the material constants cannot be 
assumed as being negligible. The foundation and contact problems in soil 
mechanics and the wave propagation problems in the earth's crust may be 
mentioned as some of the examples. The second group of problems relates to 
the fracture of essentially nonhomogeneous solids. "Hydraulic fracturing" 
of the medium which consists of sandstone and shale, the fracture of struc- 
tural materials with periodically varying material properties (as in certain 
laminated structures), and the fracture of variety of fuse-bonded materials 
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used in electronics industry are some typical examples. The distinguishing 
feature of these materials is that the material constants are continuous and 
generally differentiable functions of the space coordinates, whereas in the 
standard particulate, layered, and fiber-reinforced composites the material 
constants are discontinuous functions. The consequence of the discontinuous 
behavior of the material in crack problems is that the nature of singularity 
of the stresses at the crack tip which is on the interface and at the point 
of intersection of a crack and the interface is quite different than the 
singularity exhibited by a crack tip which is fully imbedded in a homogeneous 
medium. Even though no systematic study of the problem appears to have been 
made, it is reasonable to expect that in nonhomogeneous materials with con- 
tinuous and continuously differentiable elastic constants the nature of the 
stress singularity at a crack tip would be identical to that of a homogeneous 
solid. The existing solutions of crack and punch problems in certain specific 
nonhomogeneous materials seem to support this view. 

In most of the existing solutions of problems relating to nonhomogeneous 
solids it is assumed that the material is isotropic, the Poisson's ratio is 
constant, and the Young's (or the shear) modulus is either an exponential or 
a power function of a space variable [1], In the wedge problem described 
in [2] the class of functions E(r,e) for the Young's modulus leading to a 
feasible solution has been investigated. Some sample studies of the 
Boussinesq and contact problems for a nonhomogeneous half space may be found 
in [3]-[8]. The corresponding "torsion" problem for a half space is described 
in [9] and [10]. The equivalent crack problems in a nonhomogeneous medium 
under torsion and under anti-plane shear loading are discussed in [11] and 
[12], respectively. In the studies described in [3] and [8], the shear 
modulus is assumed to be y^expCvy) and in [12] PQexp(x“+y^) , where y=0 
is either the boundary of the half plane or the plane of the crack. In 
[4]-[7] and [9]-[ll] it is assumed that y = PqIyT* (0<m<l). This latter 
assumption clearly has the undesirable physical feature in that at y=0 the 
shear modulus becomes zero. It is particularly difficult to attach any 
physical meaning to the solution of a crack problem carried out under this 
assumption. This difficulty has been removed in the crack problem considered 
in [13] where it was assumed that the shear modulus is given by y=yQ/(l+c[y | ) , 
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where c is a constant. It should again be noted that the solutions given 
in [3]-[13] are based on the assumption that the Poisson's ratio is constant. 

In the studies mentioned above (with the exception of [12] which deals 
with the relatively simple problem of anti-plane shear) it is assumed that 
in the direction(s) parallel to the boundary of the half plane or the plane 
of the crack the shear modulus does not vary. In the case of fracture of 
such a nonhomogeneous medium since, generally, the plane of the crack is 
not a plane of symmetry, the in-plane shear component of the stress inten- 
sity factor would not be zero and hence the propagating crack would even- 
tually align itself parallel to the direction in which the modulus varies. 

In the problem considered in this paper it is then assumed that the crack 
is located on the y=0 plane, the Young's modulus is an exponential function 
of X, and the Poisson's ratio is constant (Fig. 1). 


2, Formulation of the Crack Problem 

Consider the plane elasticity problem for a nonhomogeneous solid in 
which the Poisson's ratio v is constant and the Young's modulus E is a 
function of x and y. Let F(x,y) be the Airy stress function. The stresses 
are given by 

- - - JiL m 

'^xx “ 3^ ’ *^yy 9x^ ’ ^xy 3xay * ^ ^ 

Substituting from (1) through the Hooke's Law into the compatibility equa- 
tion, for the plane problem we obtain 

- 2E(f ^ * f |:)V2F t 2(Hv)(2 || f - E ^ 

" 0 " '■E 0 ] If = “ • 


- 3 - 



Equation (2) Is for the generalized plane stress. The differential equation 
for plane strain is obtained by replacing E and v by E/(l-v2) and v/(l-v), 
respectively. From (2) it may easily be verified that if we let 


E = Eq exp(Bx+Yy) , v = constant (3) 

the differential equation becomes one of constant coefficients which may be 
written as 


V“F - 2(6 ^ + Y ^)v2F + (62-yy2) 0 + 2{1+v)6y ^ 

+ (y^-v 62) Ip- = 0 . (4) 

In problems involving the study of localized phenomena such as pertur- 
bation in stress state due to the presence of a crack or a punch a material 
representation such as (3) would not be very unrealistic. In most cases a 
reasonable approximation to the actual distribution of E(x,y) can be obtained 
by adjusting the constants E^, 3* and y. Referring to Fig, 1, in this prob- 
lem we will further assume that E is independent of y. Thus, y=0 is a plane 
of symmetry provided we also consider only those external loads which are 
symmetric with respect to y=0. It is therefore sufficient to consider one 
half of the medium, (-®<x«», y>0) only. By letting I’n (4) the differ- 
ential equation of the problem becomes 


’"P - 26(0 * 1^) * 6= 



(5) 


Note that (5) reduces to the standard bi harmonic equation for 3=0. 

Assuming the solution of (5) in the form 

GO 

F(x,y) = ^ j f(y,a)e"’^“ da , (-«<x«», y>0) , (6) 

—CO 

we obtain 


+ (2i3a-2a^-3^v) (a‘*-2i3a^-3^a^)f = 0 
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( 8 ) 


If we now look for a solution of (7) of the form f = exp(my) we find 

ni^ + (2i3a-2a2-B2v)m2 + (a'*-2i6a3-B2a2) = 0 . 

The solution of (8) is found to be 

"’l “ ""’3 ■ [(-y-|+Y 2 )/ 2 ]^ » nig = = C(-Yi“Y 2 )/ 2 ]^ » 

Y-| 2iBa-2a2-32v, = (6‘*v^-4i3^va+432va2)^ , (9) 

In (9) the roots ra^ are ordered in such a way that Re(m-j)>0, Re(m2)>0. 

It is assumed that the problem in the absence of the crack has been solved 
under the actual loading conditions, and that the crack length 2a is "small" 
compared to other (planar) dimensions of the solid. Through a superposition 
the singular part of the solution may then be reduced to that of an infinite 
nonhomogeneous plane with the self-equilibrating crack surface tractions 
as the only external loads. Thus, in the problem of interest the stresses 
and displacements vanish as (x^+y^).,*, and the solution of (7) may be 
expressed as 

-m,y ^ 

f(yjci) = A-|(a)e + A 2 (a)e , (0<y«») . (10) 

From (1), (6) and (9) it then follows that 

00 

1 r 2 

—CO 

00 

<iyy(x.y) - - ^ I a2 E Aje da , (12) 

—00 

' 2 “m.y 

- ST J “ f ^ . 

—00 

This completes the formulation of the problem for the half plane y>0 in 
which the functions A-j and Ag are determined from the two boundary conditions 
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at y=0, -»<x«». For example, let the half plane be subjected to trac- 
tions 


Oyy(x,0) = o(x) , a^y(x,0) = t(x) , (-»<x«») (14) 

on the boundary y=0 and be kept in equilibrium by a resultant force applied 
to the medium at infinity which is collinear with a force defined by the 
following x and y-components: 


P 


X 


00 

j T(x)dx , 

••00 



(15) 


From (12) and (13) and Ag may then be obtained as follows: 


A-|(a) 

P^2^o.) 

where 

Q-j (a) 



( 16 ) 

(17) 

(18) 


3. The Integral Equation 

We assume that the original cracked solid is loaded symmetrically 
in such a way that 


Oxy(^>0) “ 0 » -“<x<» . 


(19) 


In the perturbation problem, in addition to (19), we then have the following 
mixed boundary condition 


Oyy(x,+0) = p(x) , (-a<x<a) , 


( 20 ) 


v(x,0) = 0 , a<|x|<« , 


( 21 ) 
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where p(x) is a known function and v is the y-component of the displace- 
ment. From (13) and (19) it follows that 


m^A-| + m2A2 = 0 . 


( 22 ) 


To obtain the second equation to determine A-j and A2 we introduce a new 
unknown function g(x) by 


g(x) = V (x,+0) , 

From (21) and (23) it is seen that g(x) = 0 for |xl>a and 


a 

I g(x)dx = 0 . 


-a 


(23) 


(24) 


(25) 


By using the Hooke's law from (11) and (12) it can be shown that 

^ v(x.y) = - j (e>2+vni,2)e 

—00 

+ (a^+vm2^)e ^]e”^^“ da , (y>0) . 

Equation (23) and (25) would then give 

■ ^ ^ (a^+vmg^)] = j g(t)e(^-''“^^ dt. (26) 

From (22) and (26) the functions A-j and A2 are now determined as follows: 

■ I ^ A3 . (27) 

-a 


By substituting from (27) and (12) into (20) we obtain the following 
integral equation to determine g(x): 
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( 28 ) 



dtf 

(e+ia)(m^^-iri2=^) 


(m^e 


•mgy 


-m.y 


-nige 


^ da = p(x) , |x| 


<a 


To investigate and to separate a possible singular part of the kernel 

in (28) the asymptotic behavior of the inner integral must be examined. 

From (9) we observe that for |a|-»« m-j-i-lal and mg-^-la]. Now, by expressing 

the inner integral in (28) as 
00 

h(x,y.t) = |K(y,a)e^'^*“^^“ da (29) 


and by noting that any singular part h may have must be due to the behavior 
of K at |a|-x», we may write h as follows: 

00 

h(x.y.t) = I [K(y,a) - Uy>a)]e’t^-=‘>“ da 


^ fK,(y.a)e' da . 
-•00 


(30) 


where K<„ is the asymptotic value of K(y,a) for large values of [al. The 
first integral in (30) is uniformly convergent and, therefore, when substitu- 
ted into (28) the limit can be put under the integral sign. It may easily 
be shown that 


K<o(y ,ct) 



(31) 


and the second integral in (28) may be expressed as 

_i u E„(t-x) 

J ^ ^ [cos(t-x)a + i sin(t-x)a]da = • 

*00 

Also, by defining 


( 32 ) 
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M(a) = 


9 


(33) 


J(t-x)a 

(g+ia)(m^+m2) 

and by substituting y=0 in the first integral, (30) may now be written as 
follows: 

h(x,y,t) = EgJ [H(a)+M(-c.) - s1n(t-x)a]da + • (34) 

0 

Finally, if we substitute from (34) into (28) and go to limit we obtain 

7 I [|t 3 ^ + k(x,t)]g(t)dt = ^ p(x) , (-a<x<a) , (35) 

-a ° 

where the Fredholm kernel is defined by 

00 

k(x,t) = e^^ j [M(a) + M(-a) - sin(t-x)a]da , (36) 

0 

and Eq/2 is replaced by 4yQ/(l+K) in order to cover both generalized plane 
stress and plane strain problems. Here is the shear modulus at x=0, 
i.e., Uq = Eq/2(1+v), k = 3-4v for plane strain and k = (3-v)/(l+v) for 
the generalized plane stress. Note that for 3 = 0, K = Koo, k(x,t) = 0, 
and (35) reduces to the known integral equation of the simple (Mode I) 
crack problem for a homogeneous plane. 

For numerical solution the interval (-a, a) is normalized by defining 

s = t/a, r = x/a, <|)(s) = g(t), n(r,s) = k(x,t), q(r) = p(x), 

-l<(r,s)<l, -a<(x,t)<a . (37) 

In terms of the normalized quantities the integral equation (35) and the 
single-valuedness condition (24) may be expressed as 
1 

1 f 14. 

7j n(r,s)](})(s)ds =^q(r) , (-l<r<l) , (38) 
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( 39 ) 


1 

j (|)(s)ds = 0 . 
-1 


4. Stress Intensity Factors 

The index of the singular integral equation (38) is +1. Therefore its 
solution is of the following form: 

= G(s]_ ^ ^ 

where G(s) is a bounded function. The unknown function G may be determined 
from (38) and (39) to any desired degree of accuracy by using a Gaussian 
integration technique to solve the singular integral equation (see, for 
example, [14]). By observing that the left-hand side of (35) gives Oyy(x,0) 
for Ix|i>a as v/ell as for lxl<a, through a simple asymptotic analysis, the Mode 
I stress intensity factors at the crack tips defined by 

k,(a) = lim /2(x-a) a (x,0) , (41) 

x-va 

k,(-a) = lim /2(-x-a) c (x,0) , (42) 

x->-a 

, . f 

may be expressed in terms of G(s) as follows^ 


k-|(a) = - G(1)»'T , (43) 

kl(-a) = • <'*'*> 

From (23) and (37) it is seen that after obtaining G(s) the crack surface 
displacement may be calculated as 


' Note that VQexp(3a) = v(a) and the expressions (43) and (44) are iden- 
tical to those found for the homogeneous materials. 
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x/a 

viya f G(sL e-sss (45) 

It should again be emphasized that the structure of the integral equa- 
tion (35) is essentially the same as that of a homogeneous medium, namely 
its kernel has a simple Cauchy singularity. Therefore, its solution and 
consequently the stress state around the crack tip would have the conventional 
square root singularity (see (40)-(42)). 

5. Results and Discussion 

The crack problem is solved for two types of loading. In the first 
it is assumed that the plane is loaded by prescribing the displacements in 
such a way that in the uncracked medium we have 

£yy(x,0) = Eg + e^(x/a) , ejjy(x,0) = 0, ° 

From (46) the crack surface tractions in the perturbation problem may be 
expressed as follows; 

■Jyy(x.O) = P(x) = ' “■ 

In the second type of loading crack surface traction p(x) will simply be 
assumed to be a polynomial of the form 

p(x) = -P„ - P,(f) - PjCf)'^ - PjCf)’ . (43) 

In the normalized integral equation (38) the material parameter B 
enters into the kernel through ag only. Thus the calculated stress inten- 
sity factors are given with as as the variable. The results are given in 
Tables 1 and 2. Note that in the nonhomogeneous material problem since 
the kernel is also dependent on the Poisson's ratio v, the solution must 
be obtained for a given value of v, and for the cases of plane stress 
and plane strain separately. Tables 1 and 2 give the normalized stress 
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intensity factors at the crack tips a and -a for plane stress and plane 
strain cases, respectively, by assuming that v = 0,3. The results are 
obtained by taking only one of the six input parameters e^, e-j, p^, p^, 

P 2 , and Pj nonzero at a time. Since the problem is linear the results can 
be superimposed in any suitable manner. Note that for 3a-^, that is, for 
a given a and 3->0 or a given 0 and a->0, as expected, the stress intensity 
factor ratios reduce to those of a homogeneous medium, which, for an arbi- 
trary traction p(x), are given by 

a 

ki(a) =-Lf p(x) . (49) 



(50) 


For 6a = 0.5 the effect of Poisson's ratio on the stress intensity 
factors in plane stress and plane strain cases is shown in Table 3. The 
results show that this effect is rather insignificant. Consequently, as 
seen from Tables 1-3 the difference between plane stress and plane strain 
results is also insignificant. 

Since the effect of the Poisson's ratio and the thickness constraint 
on the stress intensity factors is negligibly small, from the results given 
in Tables 1 and 2 it is possible to develop simple empirical formulas for 
the stress intensity factors. For example, from Fig, 2 reproducing the 
results for a uniformly pressurized crack (columns 6 and 7, Table 1) it may 
be seen that the stress intensity factors vary approximately linearly with 
the variable ea. Hence, in this case, the following approximate formulas 
may be used to evaluate the stress intensity factors; 


(51) 


ki(a) = Pq*^( 1 + 0.22ea) , 
k-|(-a) = Pq*^(^ " 0.266a) . 
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A sample result showing the crack surface displacement y(x) for a uni- 
formly pressurized crack as calculated from (45) is given in Fig, 3. The 
figure also shows v(x) for a homogeneous material. As expected, in the 
nonhomogeneous medium in the stiffer portion of the material the crack 
surface displacement is smaller than that of the homogeneous medium, and 
the reverse trend may be observed in the less stiff portion of the material. 
In this as well as in the previous studies the assumption of constant 
Poisson's ratio has been made for analytical reasons. Even though not very 
conclusive, the results given in this paper show that neglecting the possi- 
ble special variation of the Poisson's ratio is not a very restrictive 
assumption. 
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Table 1. The normalized stress intensity factors for various loading 
conditions for the case of generalized plane stress (v=0.3). 


3a 

EHfOMHfgll 



k^(-a) 

k^(a) 

Hlgll 

HiBMiaieai 

IHtQMiaigll 


®o^o*^ 


eiEo*^ 

Po*^ 

Po'^ 


Pi»^ 

Pj^ 

Pg*^ 

Ps*^ 

Pj^- 

-0 

1.0 

1.0 

0.5 

-0.5 

1.0 

1.0 

0.5 

-0.5 

0.5 

0.5 

0.375 

-0.375 

0.01 

1.008 

0.992 

0.505 

-0.495 

1.003 

0.997 

0.500 

-0.500 

0.501 

0.499 

0.375 

-0.375 

0.10 

1.078 

0.925 

0.552 

-0.453 

1.025 

0.973 

0.500 

-0.500 

0.506 

0.493 

0.375 

-0.375 

0.25 

1.202 

0.820 

0.640 

-0.389 

1.060 

0.930 

0.498 

-0.499 

0.515 

0.483 

0.374 

-0.375 

0.50 

1.435 

0.665 

0.814 

-0.302 

1.113 

0.861 

0.495 

-0.495 

0.528 

0.466 

0.372 

-0.373 

0.75 

1.713 

0.535 

1.031 

-0.234 

1.162 

0.797 

0.489 

-0.489 

0.540 

0.450 

0.369 

-0.370 

1.00 

2.048 

0.429 

1.301 

-0.181 

1.209 

0.740 

0.483 

-0.480 

0.552 

0.435 

0.366 

-0.366 
























Table 2. The normalized stress intensity factors for various loading conditions 
for the case of plane strain (v=0.3). 


3a 

n-v2)k,(a) 

(l-v2)k^(-a) 

(l-v2)k^(a) 

0-v^k,(-a) 

k^(a) 

k^(-a) 

k^(a) 

k^(-a) 

k^(a) 

k^(-a) 

k^(a) 

k-|(-a) 

^0^0*^ 



^1^0^ 

Po>^ 

Po-^ 



P2*^ 

P2*^ 

Pj^- 

P3^ 

-*•0 

1.0 

1.0 

0.5 

-0.5 

1.0 

1.0 

0.5 

-0.5 

0.5 

0.5 

0.375 

-0.375 

0.01 

1.008 

0.992 

0.505 

-0.495 

1.003 

0.997 

0.500 

-0.500 

0.501 

0.499 

0.375 

-0.375 

0.10 

1.078 

0.925 

0.552 

-0.453 

1.026 

0.973 

0.500 

-0.500 

0.506 

0.493 

0.375 

-0.375 

0.25 

1.203 

0.821 

0.640 

-0.389 

1.061 

0.931 

0.498 

-0.499 

0.515 

0.483 

0.374 

-0.375 

0.50 

1.439 

0.667 

0.814 

-0.302 

1.117 

0.863 

0.494 

-0.495 

0.529 

0.466 

0.372 

-0.373 

0.75 

1.721 

0.539 

1.032 

-0.234 

1.170 

0.801 

0.489 

-0.489 

0.542 

0.451 

0.369 

-0.370 

1.00 

2.063 

0.433 

1.304 

-0.181 

1.222 

0.745 

0.483 

-0.481 

0.555 

0.436 

0.366 

-0.366 







Table 3. The effect of Poisson's ratio on the stress intensity factors, 
(ea = 0.5) (Eq* = Eq for plane stress, Eq* = 

for plane strain). 


■ 



■aigiM 

lasa 

manmw 

mai 

am 

tsmrn 

lasi 

msi 

BUBI 

IIBSI 

la&i 

■ 

V 

0 0 

0 0 



Po-^ 



Pl»^ 

P2»^ 


Ps-^ 

PjV^ 

(/) 

CO 

0.01 

1.425 

0.660 

0.813 

-0.302 

1.104 

0.856 

0.494 

-0,494 

0.525 

0.464 

0.372 

-0.372 

u 

4-> 

(/) 

0.15 

1.430 

0.662 

0.813 

-0.302 
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Fig. 1. The crack geometry in the nonhomogeneous medium and the 
variation of the Young's modulus E = 


y 



Fig.; 2. Stress intensity factors in a nonhomogeneous medium having 
a uniformly pressurized crack (E(x) = v = 0,3, plane 

stress conditions). 




x/a 


Fig. 3. Crack surface displacement v(x) in a nonhomogeneous and a 
homogeneous medium under uniform pressure p^ applied to the 
crack surfaces, (E(x) = EQexp(0,5x/a) for the nonhomogeneous 
medium, E(x) = E^ for the homogeneous medium, v = 0,5, pa » 
0.5, plane stress conditions). 
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